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r-| ! Abstract 

We present a consistent calculation of bubble-nucleation rates in theories of two 
rS \ scalar fields. Our approach is based on the notion of a coarse-grained free energy 

that incorporates the effects of fluctuations with momenta above a given scale k. 
We establish the reliability of the method for a variety of two-scalar models and 
confirm the conclusions of previous studies in one-field theories: Langer's theory of 
homogeneous nucleation is applicable as long as the expansion around the semiclas- 
sical saddle point associated with tunnelling is convergent. This expansion breaks 
down when the exponential suppression of the rate by the saddle-point action be- 
comes comparable to the pre-exponential factor associated with fluctuations around 
the saddle point. We reconfirm that Langer's theory is not applicable to the case 
of weakly first-oder phase transitions. We also find that the same is true in gen- 
eral for radiatively induced first-order phase transitions. We discuss the relevance of 
our results for the electroweak phase transition and the metastability bound on the 
Higgs-boson mass. 



1 Introduction 



The calculation of bubble-nucleation rates during first-order phase transitions is a difficult problem, 
both at the conceptual and at the technical level. The standard approach to this problem is based 
on the work of Langer 0. His formalism has been applied to relativistic field theory by Coleman Q 
and Callan || and extended by Affleck Q and Linde Q to finite temperature. The nucleation rate I 
gives the probability per unit time and volume to nucleate a certain region of the stable phase (the 
true vacuum) within the metastable phase (the false vacuum). Its calculation relies on a semiclassical 
approximation around a dominant saddle-point, which is identified with the critical bubble. This is a 
static configuration (usually assumed to be spherically symmetric) within the metastable phase whose 
interior consists of the stable phase. It has a certain radius that can be determined from the parameters 
of the underlying theory. Bubbles slightly larger than the critical one expand rapidly, thus converting 
the metastable phase into the stable one. The nucleation rate is exponentially suppressed by the 
action (the free energy rescaled by the temperature) of the critical bubble. Possible deformations of 
the critical bubble generate a pre-exponential factor. The leading contribution to this factor has the 
form of a fluctuation determinant and corresponds to the first-order correction to the semiclassical 
result. 

For a four-dimensional theory of a real scalar field at temperature T, in the limit that thermal 
fluctuations dominate over quantum fluctuations, the bubble-nucleation rate is given by f|]-P] 
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Here V is the free energy of the system for a given configuration of the field (ft that acts as the 
order parameter of the problem. The rescaled free energy of the critical bubble is S = Tb/T = 
[r (^>b(0) ~~ -T(O)] /T, where 0b( r ) is the spherically-symmetric bubble configuration and = corre- 
sponds to the false vacuum. The fluctuation determinants are evaluated either at at = or around 
4> = 0b( r )- The prime in the fluctuation determinant around the bubble denotes that the three zero 
eigenvalues of the operator [5 2 T / Scj) 2 ]^^ have been removed. Their contribution generates the fac- 
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tor (S/2tt) ' and the volume factor that is absorbed in the definition of I (nucleation rate per unit 
volume). The quantity Eq is the square root of the absolute value of the unique negative eigenvalue. 

In field theory, the free energy (more precisely the thermodynamic potential) density T of a system 
for homogeneous configurations is usually identified with the temperature-dependent effective poten- 
tial. This is evaluated through some perturbative scheme, such as the loop expansion Q. The profile 
and the free energy of the bubble are determined through the potential. This approach, however, faces 
three fundamental difficulties: 

a) The effective potential, being the Legendre transform of the generating functional for the con- 
nected Green functions, is a convex function of the field. Consequently, it does not seem to be 
the appropriate quantity for the study of tunnelling, as no structure with more than one minima 
separated by a barrier exists. It has been argued in ref. J/J that the appropriate quantity for the 
study of tunnelling is the generating functional of the 1PI Green functions (calculated pertur- 
batively), which differs from the effective potential in the non-convex regions. However, as we 
discuss in the following, the consistent picture must rely on the notion of coarse graining and on 
the separation of the high-frequency fluctuations that may be responsible for the non-convexity 
of the potential, from the low-frequency ones that are relevant for tunnelling. Such notions 
cannot be easily implemented in the context of perturbation theory. 

b) The fluctuation determinants in the expression for the nucleation rate have a form completely 
analogous to the one-loop correction to the potential. The question of double-counting the effect 
of fluctuations (in the potential and the prefactor) must be properly addressed. This point is 
particularly important in the case of radiatively induced first-order phase transitions. These 
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are a consequence of the appearance of a new vacuum state in the theory as a result of the 
integration of (quantum or thermal) fluctuations H. A radiatively induced first-order phase 
transition takes place in theories for which the tree-level potential has only one minimum, while 
a second minimum appears at the level of radiative corrections^]. 

Another difficulty concerns the ultraviolet divergences that are inherent in the calculation of 
the fluctuation determinants in the prefactor. An appropriate regularization scheme must be 
employed in order to control them [j^|- [jl2^ . Moreover, this scheme must be consistent with the 
one employed for the absorption of the divergences appearing in the calculation of the potential 
that determines the free energy of the critical bubble. 



In a previous publication [|13|], based on the framework described in refs. [14, [Tq ], we demonstrated 
that all the above issues can be resolved through the implemention of the notion of coarse graining in 
the formalism. As the appropriate quantity for the description of the physical system we employed the 
effective average action |lq| , which is the generalization in the continuum of the blockspin action 
of Kadanoff []17|. It can be interpreted as a coarse-grained free energy at a given scale k. Fluctuations 
with characteristic momenta q 2 i; k 2 are integrated out and their effect is incorporated in IV In 
the limit k — > 0, T/% becomes equal to the effective action. The k dependence of is described by 



an exact flow equation [18], typical of the Wilson approach to the renormalization group [19]. This 



flow equation can be translated into evolution equations for the invariants appearing in a derivative 



expansion of the action [20, 21 1. In ref. [13] we considered only the effective average potential Uk and a 
standard kinetic term and neglected higher derivative terms in the action. We shall employ the same 
approximation in this paper also. Its validity is guaranteed by the small anomalous dimensions in the 
models we consider. The bare theory is defined at some high scale A that can be identified with the 
ultraviolet cutoff. At scales k below the temperature T, the theory can be described in terms of an 
effective three-dimensional action at zero temperature |23|]. This dimensional reduction indicates 
the absence of explicit time dependence for the parameters of the theory at low energy scales. 



In ref. [13] we considered as a starting point the action Tk for a real scalar field at a scale 
ko below the temperature, such that the theory has an effective three-dimensional description. We 
approximated by a standard kinetic term and a potential with two minima. We assumed that 
this form of the potential results from the bare potential U\ after the integration of (quantum and 
thermal) fluctuations between the scales A and ko. Some of these fluctuations may correspond to 
additional massive degrees of freedom that decoupled above the scale ko. We computed the form of 
the potential at scales k < ko by integrating an evolution equation derived from the exact flow 
equation for IV is non-convex for non-zero k, and approaches convexity only in the limit k — > 0. 
The nucleation rate must be computed for k larger than the scale kf at which the functional integral 
in the definition of starts receiving contributions from field configurations that interpolate between 
the two minima. This happens when — k 2 becomes approximately equal to the negative curvature at 
the top of the barrier |24j|. For k > kf the typical length scale of a thick- wall critical bubble is ^ 1/k. 
Through the use of Uk , the first problem in the calculation of the nucleation rate mentioned above is 
resolved. 

We performed the calculation of the nucleation rate for a range of scales above and near kf. In 
our approach the pre-exponential factor is well-defined and finite, as an ultraviolet cutoff of order k is 
implemented in the calculation of the fluctuation determinants, such that fluctuations with character- 
istic momenta q 2 ^ k 2 are not included. This is a natural consequence of the fact that all fluctuations 



1 In ref. |8| an alternative procedure was suggested for the treatment of radiatively-induced first-order phase transi- 
tions: The fields whose fluctuations are responsible for the appearance of the new vacuum are integrated out first, so 
that an "effective" potential with two minima is generated for the remaining fields. Our philosophy is different: We 
integrate out high-frequency fluctuations of all fields, so that we obtain an effective low-energy action which we use 
for the calculation of the nucleation rate. Our procedure involves an explicit infrared cutoff in the calculation of the 
low-energy action. This prevents the appearance of non-localities arising from integrating out massless fields, which may 
be problematic for the approach of ref. bl . For example, the fields that generate the new vacuum in radiatively-induced 
first-order phase transitions are usually massless or very light at the origin of the potential. 
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with typical momenta above k are already incorporated in the form of U}.. This modification also 
resolves naturally the problem of double-counting the effect of the fluctuations. 

However, an important issue arises at this point. The scale k was introduced in the problem 
as a mere calculational tool. If our approach makes sense, the choice of k should not affect physical 
parameters such as the nucleation rate. The remarkable outcome of our study was that this expectation 
was confirmed. We found that the saddle-point configuration has an action Sk with a significant 
k dependence. For strongly first-order phase transitions, the nucleation rate I = Afcexp(— Sk) is 
dominated by the exponential suppression. The main role of the prefactor Ak, which is also k- 
dependent, is to remove the scale dependence from the total nucleation rate. The implication of our 
results is that the critical bubble should not be identified just with the saddle point of the semiclassical 
approximation. It is the combination of the saddle point and its possible deformations in the thermal 
bath (accounted for by the fluctuation determinant in the prefactor) that has physical meaning. We 
also found that, for progressively more weakly first-order phase transitions, the difference between 
Sk and ln(Afc//c|) diminishes. This indicates that the effects of fluctuations become more and more 
enhanced. At the same time, a significant k dependence of the predicted nucleation rate develops. 
The reason for the above deficiency is clear. When the nucleation rate is roughly equal to or smaller 
than the contribution from the prefactor, the effect of the next order in the expansion around the 
saddle point is important and can no longer be neglected. This indicates that there is a limit for the 
validity of Langer's picture of homogeneous nucleation |Q] . 

In this work we describe how our method can be applied to a more complicated system. As such 
we have chosen a theory of two scalar fields. It provides a framework within which we can test the 
reliability of our computation of the nucleation rate in the case of two fluctuating fields. The evolution 
equation for the potential resembles very closely the ones appearing in gauged Higgs theories, with 
the additional advantage that the approximations needed in the derivation of this equation are more 
transparent. Moreover, we expect the qualitative conclusions for the region of validity of Langer's 
picture of homogeneous nucleation to be valid also for gauged Higgs theories. The most interesting 
feature of the two-scalar models is the presence of radiatively induced first-order phase transitions. 
Such transitions usually take place when the mass of a certain field is generated through the expectation 
value of another. The fluctuations of the first field can induce the appearance of new minima in the 
potential of the second, resulting in first-order phase transitions [0] . As we have already discussed, 
the problem of double-counting the effect of fluctuations is particularly acute in such situations. 
The introduction of a coarse-graining scale k resolves this problem, by separating the high-frequency 
fluctuations of the system which may be responsible for the presence of the second minimum through 
the Coleman- Weinberg mechanism, from the low-frequency ones which are relevant for tunnelling. 

In the following sections we present the calculation of the bubble-nucleation rate for first-order 
phase transitions in a theory of two real scalar fields. We do not discuss the evolution of for 
k^T. We start the evolution at a scale k$ sufficiently below the temperature of the system, so that 
the dynamics is three-dimensional to a good approximation^]. As an initial condition we consider a 
potential Uk , whose form is determined by the bare potential U\ and the integration of fluctuations 
between the scales A and fco. We first establish the reliability of our approach by considering potentials 
Uk with two minima. Then we turn to the radiatively induced first-order phase transitions, for which 
the second minimum is generated at some point in the evolution at a scale k < ko. In all cases, we 
integrate the evolution equation for the effective three-dimensional theory starting at the scale /co, and 
perform the calculation of the nucleation rate as described earlier in the introduction. 

In the following section we present the model we consider and derive the evolution equation for the 
potential. In section 3 we summarize the technical points in the calculation of the nucleation rate. Our 
results are presented in section 4. The implications of our results for the first-order phase transitions 
in gauged Higgs models (such as the electroweak phase transition) are given in the conclusions of 

2 For readers who are interested in the details of the mechanism of dimensional reduction in our approach, detailed 
discussions can be found in refs. pi B3L E6[ for a variety of models. 
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section 5. 



2 The model and the evolution equation for the potential 



We consider a model of two real scalar fields <f>\ and 4>2- The effective average action Tk((pi,(j) 2 ) [|16| 
results from the effective integration of degrees of freedom with characteristic momenta larger than 
a given cutoff k. This is achieved by adding an infrared cutoff term to the bare action, so that the 
effective action does not receive contributions from modes with characteristic momenta q 2 ^ k 2 . In 
this work we use the simplest choice of a mass-like cutoff term ~ k 2 (cj) 2 +</>2), for which the perturbative 
inverse propagator for massless fields is Pk(q) ~ q 2 +k 2 . Subsequently, the generating functional for the 
connected Green functions is defined, from which the generating functional for the 1PI Green functions 
can be obtained through a Legendre transformation. The presence of the modified propagator in the 
above definitions results in the effective integration of only the fluctuations with q 2 ^ k 2 . Finally, the 
effective average action is obtained by removing the infrared cutoff from the generating functional for 
the 1PI Green functions. 

The effective average action obeys an exact flow equation, which describes its response to 
variations of the infrared cutoff k |l8| . This can be turned into evolution equations for the invariants 



appearing in a derivative expansion of [20]. In this work we use an approximation which neglects 
higher derivative terms in the action and approximates it by 



(2.1) 



The above action describes the effective three-dimensional theory that results from the dimensional 
reduction of a high-temperature four-dimensional theory at scales below the temperature. The tem- 
perature has been absorbed in a redefinition of the fields and their potential, so that these have 
dimensions appropriate for an effective three-dimensional theory. The correspondence between the 
quantities we use and the ones of the four-dimensional theory is given by 



'1,24 



; 1,2 



U((f>U<t>2 



(2.2) 



In this way, the temperature does not appear explicitly in our expressions. This has the additional 
advantage of permitting the straightforward application of our results to the problem of quantum 



tunnelling in a three-dimensional theory at zero temperature [13]. 



The evolution equation for the potential can be written in the form [25 



1 

8^ 



k 2 + M?(<k,<h) ~ \/k 2 + M 2 (0, 0) + 
+ ^k 2 + M|(0!,0 2 ) - \A 2 + M 2 2 (0,0) 



(2.3) 



where M 2 2 {4>i, $2) are the two eigenvalues of the field-dependent mass matrix, given by 



(U u - U22? + 4^ 



(2.4) 



%J = d 2 U k / d(j)id(j)j. The only neglected corrections to eq. (2^3) are related to the wave-function 



with U, 

renormalization of the fields. We expect these correction to be small, as the anomalous dimension 
is r] ~ 0.035 — 0.04 for the models we consider and the evolution of the potential takes place over a 
limited range of k. We consider models with the symmetry (f>2 «-> —02 throughout this paper. This 
means that the expressions for the mass eigenvalues simplify along the 1^1-axis: M 2 = d 2 Uk/d(j) 2 , 
M| = d 2 Uk /d(j) 2 ■ We point out that our method does not face the problems of non-convergence 
of perturbation theory that often appear in the context of radiative symmetry breaking ||. Up to 
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wave-function renormalization effects, eq. ( |2.3| ) is exact and its numerical solution is very accurate (see 
below). However, the inclusion of higher-derivative terms in the action is required for the quantitative 
study of models with strongly coupled phases, such as the symmetric phase of the electroweak theory 
above the critical temperature of the electroweak phase transition. 

The first step of an iterative solution of eq. (^3) gives [27] 



+ 



+ 



U ko (fa,fa)-U ko (0) + 

det[-d 2 + k 2 + M 2 (fa,fa)] det[-8 2 + k 2 + M 2 (0)] 
det[-<9 2 + k 2 + Mf(fa, fa)] det[-d 2 + k 2 + M 2 (0)] 

det[-8 2 + k 2 + M 2 (fa,fa)] det[-d 2 + k 2 + Mf (0)] 
det[-d 2 + k 2 + M 2 (fa, fa)] det[-d 2 + k 2 + M|(0)] 



2 



1 



In 



(2.5) 



For k — > this is a regularized one-loop approximation to the effective potential. Due to the ratio of 
determinants, only momentum modes with k 2 ^ q 2 Si fcg are effectively included in the momentum 
integrals in eq. ( |2.5| ). The above expression demonstrates the form of ultraviolet regularization of 
fluctuation determinants that is consistent with the cutoff procedure that we described in the beginning 
of this section in the derivation of the evolution equation for the potential. An analogous regularization 
will be used in the following section for the fluctuation determinants in the expression for the nucleation 



rate. The solution of eq. (2J3) also demonstrates that the decoupling of heavy modes is automatically 
built in the evolution equation for the potential. If M 2 (fa,fa) 3> k^ > k 2 , the contribution of the 
respective mode to this solution is negligible. 

We can determine the symmetries of the theory by specifying the form of the potential at the scale 
ko < T at which we start the evolution. As we have already mentioned, we assume the symmetry fa <-> 
—fa throughout this paper. In section 4 we consider an example of a theory without any additional 
symmetry for the fa field. However, for most of our results and especially for those concerning 
radiatively induced first-order phase transitions, we use potentials invariant under fa <-> — fa, fa <-> 
—fa, fa <-> fa. The solutions of eq. ( |2.3| ) preserve the above symmetries for k < ko- 

For the numerical integration of eq. fl2.3| ) 



we use the algorithms described in ref. [28|, with an 
appropriate generalization to the two-field case. As we shall be interested in the form of the potential 
along the fa-axis, we perform the integration of the evolution equation for Uk(fa, fa) in a thin region 
around this axis. In certain cases, which would require excessive computer time, we solve approximate 
evolution equations for Uk and d 2 Uk/d(j) 2 . along the 0i-axis (fa = 0). We control the accuracy of these 
approximate solutions by comparing them against the full solution for U)-(fa, fa) when this is possible. 
Moreover, we check our numerical results against exact analytical solutions for simplified forms of the 
evolution equation pS[ |. We estimate that the numerical uncertainties introduce an error of order 1% 
to our results, smaller than the error generated by the omission of the anomalous dimension. 



3 The calculation of the bubble-nucleation rate 

The calculation of the nucleation rate proceeds in complete analogy to the one described in detail 
in ref. |l3| for the one-field case. For this reason, we summarize here only the main steps of the 
calculation. We consider potentials U k (fa, fa) with two minima along the c^i-axis (fa = 0): the stable 
(true) minimum is located at fa = fa and the unstable (false) one at fa = <f>f = 0. 

The nucleation rate is exponentially suppressed by the action Sk (the rescaled free energy) of the 
saddle-point configuration fa(r) that is associated with tunnelling. This is an SO(3)-invariant solution 
of the classical equations of motion along the 0i-axis (fa = 0) which interpolates between the local 
maxima of the potential —Uk(fa) = —Uk(fa, fa = 0). It satisfies the equation 

d 2 fa | 2 dfa _ dU k (fa) 
dr 2 r dr dfa ' 
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with the boundary conditions < 
saddle point is given by 

S k = 4vr 



!>b — ► for r — > oo and d(f>\,/dr = for r = 0. The action of the 



dr 



+ £4(0b(r)) - C/ fc (0) 



r 2 dr. 



(3.2) 



The profile of the saddle point can be easily computed with the "shooting" method |3C]. We integrate 
eq. (|3.1|) numerically, starting at r = with a value of 0i near the true minimum 0< and d<j)\/dr = 0. 
We then adjust the initial value of 0i so that the boundary condition 0t, — > for r — > oo is satisfied. 

In all the models we consider, the minima of the potential and the saddle-point configuration are 
located along the 0i-axis (02 = 0). The unrenormalized decay rate per unit volume from the false 
minimum towards the true one is then given by [0, 0, 0] 



2tt V2vr7 



det^ry^^ det^ry^ 2 ]^^ 



det[5 2 T k /5<pl 



01=0 



det[5 2 iy<50 2 



2J0l=O 



-1/2 



exp(-Sfc). 



(3.3) 



This is analogous to eq. ([O]) after the absorption of the explicit factors of T in the redefinition of 
the fields and potential and the introduction of a coarse-graining scale. If there are several equivalent 
true vacua, the above rate must be multiplied by an appropriate factor, in order to take into account 
the possibility of the false vacuum decaying into any of them. 

The pre-exponential factor corresponds to the first correction to the semiclassical approximation 
in the saddle-point method. We are considering models with the symmetry 02 «-> —02 throughout 
this paper. This guarantees that the saddle-point configuration is located along the 0i-axis, where 
the mass eigenvalues simplify: M 2 = d 2 Uk/dcf) 2 = U\\, M| = d 2 Uk/dcf) 2 = 1/22- The numerators in 
eq. (|3.3| ) are the fluctuation determinants around the saddle-point 



det' 
det 



5 2 T k /54>\ 
5 2 T k /5<i?2 



det' 
det 



-d 2 + U ll {<p l =0 b (r)) 
-d 2 + U 2 2^i = 0b(O) 



(3.4) 



while the denominators are the fluctuation determinants around the false vacuum 



det 
det 



5 2 T k 



5 2 T> 



</>i=0 
01=0 



det 
det 



-3 2 + [/n(0i = 0) 
-3 2 + [/ 22 (0i = 0) 



(3.5) 



The differential operator —d 2 + Un(<pb{ r )) has three zero modes (the three spatial translations of the 
critical bubble). The prime over the determinant indicates that these modes have to be omitted in its 
calculation. Their contribution generates the factor (Sfc/27r) 3//2 in eq. ( |3,3j ) and the volume factor that 
is absorbed in the definition of I (nucleation rate per unit volume). The quantity Eq is the square 
root of the absolute value of the unique negative eigenvalue of the above operator. The operator 
— d 2 + U22(4>h{ r )) has only positive eigenvalues, because the bubble is stable in the 02 direction at 
02 = 0. 

The pre-exponential factor defined in eq. (^^) is in general ultraviolet-divergent and an appropriate 
regularization scheme must be employed. Within our approach, the form of the regularization is 
dictated by the discussion at the end of the previous section. The effect of the high-frequency modes 
has been incorporated in the form of the coarse-grained potential Uk, which is obtained through the 
integration of the evolution equation ■ Fluctuation determinants computed within the low-energy 
theory must be replaced by a ratio of determinants, in complete analogy to eq. (p.5|). This implies 
that the nucleation rate is given by 



^i fc ^2fcexp(-5 fc ) 



where 
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Eo (Sk\ 3/2 
2vr \2n) 



det' [-a 2 + U u (Mr))] det [-d 2 + /c 2 + C/n 



det [-d 2 + k 2 + Unfair))] det [-5 2 + U n (0)] 



-1/2 
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2 A: 



det [-a 2 + U 22 fa(r))] det [-<9 2 + A; 2 + t/ 22 (0)] 



det [-a 2 + k 2 + U 22 fa(r))} det [-8 2 + f/ 22 (0)] 



-1/2 



(3.6) 



The above form of the pre-exponential factors guarantees that only modes with characteristic momenta 
q 2 ^ k 2 contribute to the nucleation rate. Another feature of eq. (|3,6j ) is the decoupling of heavy modes 
with M 2 > k 2 . 

The differential operators that appear in eq. (|3.6|) have the general form 



-d 2 + m 2 K + aW ik (r) 



where 



Uii(0)+Kk 2 , 

W lk (r) = Unfair)) - U ti (0), 



(3.7a) 

(3.7b) 
(3.7c) 



with i = 1 or 2 and k, a = or 1. Since the Wj KCt operators are SO(3) symmetric, it is convenient to 
use spherical coordinates and express the eigenfunctions in terms of spherical harmonics. This leads 
to 



det W iKa 



\2m 



ntdetnW*) 5 
-V 2 e +m 2 K + aW ik (r), 



where 



V 2 = 



dr 2 



£(£ + 1) 



(3.8) 
(3.9) 



and £ is the usual angular quantum number. 

The computation of such complicated determinants is made possible by a theorem [31, [| that 
relates ratios of determinants to solutions of ordinary differential equations. In particular, we have 



9Uk 



det HW _ det[-V| + mf K + 1 • W ik (r)] _ yu Kl (r -> cxd) 



det Wu K o det[- Vf + m 2 lK + • W ik {r)] y ieK0 {r 
where yunair) is the solution of the differential equation 



oo 



(3.10) 



d2 W + l ) 2 
-^2+^2-^+ m l + aW *k(r) 



yuna(r) 



(3.11) 



with the behaviour yu Ka {r) oc r e+l for r — > 0. Such equations can be easily solved numerically with 
Mathematica [32]. The final expression for the nucleation rate, appropriate for an efficient numerical 
computation |i"3[j , is 



JL ^V /2 
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(3.12) 



7 



(0a) C4(<fc, 0) 




1 2 3 4 5 




1 2 3 4 5 



(0c) Wifc(r 




2 4 6 
(|d) W 2fc (7 




2 4 



8 10 




2 4 6 8 10 
(@f) Results 

* * * I 
1 1 * 



50 


-50 
-100 
-150 ♦ ♦ 

-200 



0.6 0.7 0.8 0.9 1 



Figure 1: The steps in the computation of the nucleation rate for a model with initial potential given 
by eq. with m\ kQ = —m\ kn = 0.1 /c 2 , Afc = gk = 0.1 ko, Jk = 0.6 k^ 2 . The calculation is 



performed between the scales ki = e 0,8 ko and kf = e 1,2 k^. All dimensionful quantities are given in 
units ofkf. In fig. |7]f we plot the saddle-point action (diamonds) , the two pref actors \a{A\k/k^) (stars) 
andln(A 2 k) (triangles), and the nucleation rate ln(I/kj) (squares) as a function of kj y/Un(4>t, 0). 



The calculation of en is slightly complicated because of the necessity to eliminate the zero eigenvalues 
in g'ii - Also the (unique) negative eigenvalue — Eq of Wiooi must be computed for the determination 



of cio- How these steps are achieved is described in ref. [13], where it is also checked that the product 
of an infinite number of terms in eq. ( p. 12 ) is finite. This can be shown by employing first-order 
perturbation theory in Wik for large £, which gives 



c ie ^l + ^ + 0(£- 4 ) 



with 



D; 



-k 1 



r 3 W ik (r)dr. 



(3.13) 



As eq. ( 3.13 ) gives a reasonable approximation of ca even for small values of 
consider, we can derive an approximate expression for the prefactors 



ln(A ik ) w sign(A)v |A| ^r- 



in all the cases we 



(3.14) 



We have checked that the above expression gives a good analytical approximation to the numerical 
results we present in the next section. 



4 Results 

As a first example, we apply our formalism to a theory with a potential that strongly resembles the 
ones used in refs. (l3|, 25 1. The potential is consistent with the symmetry (f>2 *-> —(f>2 and near the 
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Figure 2: As in fig. \Q 6u£ /or a model with initial potential given by eq. (4^) with <f>Q k 
Afc = 0.4 ko, gk = 0.3 A;o and Vk = 1- The calculation is performed between the scales hi - 
and kj = e~°' s ko. All dimensionful quantities are given in units of kf. 



2 ko, 
-° A k 



ii-axis has the form 



1 O ,0.1 9 , 1 



U ko (^h) = -Jkji + \rn\ k ^\ + -ml ko <p 2 2 + -A feo [4>\ + 4) + 9k 4>\4>l (4.1) 

The term linear in (pi can be removed through an appropriate shift of <pi |33|, [l3[] . This would introduce 
additional terms ~ 4>\ and ~ 4>i<p2- In fig. la we present the evolution of Uk{4>i) = £4(</>i>0) for 

m iko = — 0-1 ^0) m 2k = 0-1 ^oj ^fco = 5fco = 0-1 and Jfc = 0-6 &o ■ We always shift the location of 
the false vacuum to zero. The evolution of U22{4>i) = <9 2 £4/<9</>2(</>i, 0) is displayed in fig. lb. The solid 
lines correspond to k- L /ko = e -0 ' 8 , while the line with longest dashes (that has the smallest barrier 
height) corresponds to kf /ko = e~ 12 . At the scale kf the negative curvature at the top of the barrier 
is slightly larger than —k 2 . This is the point in the evolution of the potential where configurations 
that interpolate between the minima start becoming relevant in the functional integral that defines the 



coarse-grained potential [24, 34]. For this reason, we stop the evolution at this point. The potential 
and the field have been normalized with respect to kf, so that they are of order 1. The profile of the 
critical bubble 4>b( r ) is plotted in fig. le in units of kf for the same sequence of scales. For k ~ kf the 
characteristic length scale of the bubble profile and 1/k are comparable. This is expected, because 
the form of the profile is determined by the barrier of the potential, whose curvature is ~ — k 2 at this 
point. This is an additional indication that we should not proceed to coarse-graining scales below kf. 
The quantities W\k{r) = {7ii(</>b(?"")) — Un(0) and W2k{r) = ^(^b^)) _ ^22(0) are plotted in figs, lc 
and Id respectively. 

Our results for the nucleation rate are presented in fig. If. The horizontal axis corresponds to 
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Figure 3: As in fig. \Q but for a model with initial potential given by eq. (I^.i ) with 4>Q k(j = 2 ko, 



^k = 0.35 ko, gk = 0.8 ko and = 0.8. The calculation is performed between the scales ki = e k 



and k 



I 



°' 7 ko. All dimensionful quantities are given in units of kt 



k j 'y/UuJj^t) > i- e - t ne ratio of the scale k to the square root of the positive curvature of the potential 
along the </>i-axis at the true vacuum. The latter quantity gives the mass of the field <f>\ at the 
absolute minimum. Typically, when k crosses below this mass (corresponding to the value 1 on 
the horizontal axis) the massive fluctuations of the fields start decoupling (in all the examples we 
present the mass of <p2 is of the same order or larger than that of <p± at the absolute minimum) 
and the evolution of the convex parts of the potential slows down and eventually stops. The dark 
diamonds give the negative of the action Sk of the saddle point at the scale k. We observe a strong k 
dependence of this quantity. The stars in fig. Id indicate the values of ln(Aik/kj) and the triangles 
those of ]n(A2k), where the two prefactors An-, A 2 k are defined in eqs. ( |3.6| ). Again a significant 
k dependence is observed. More specifically the values of Aik and A 2 k decrease for decreasing k. 
This is expected, because k acts as the effective ultraviolet cutoff in the calculation of the fluctuation 
determinants. For smaller k, fewer fluctuations with wavelengths above an increasing length scale 
~ contribute explicitly to the fluctuation determinants. We notice that both A^ and A 2 k are 
significantly smaller than 5^ and enhance the total nucleation rate. The dark squares give our results 
for ln(I/kf) = —Sk + ln^i^^fe/^/)- This quantity has a very small k dependence, which confirms 
our expectation that the nucleation rate should be independent of the scale k. The small residual 
dependence on k can be used to estimate the contribution of the next order in the expansion around the 
saddle point. This contribution is expected to be smaller than the first-order correction \n(AikA 2 k/kf). 
This example confirms the conclusions of refs. |l3|, on the region of validity of Langer's theory of 
homogeneous nucleation: This theory is applicable as long as the expansion around the semiclassical 
saddle-point is convergent. 
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Figure 4: ^4s in fig. 6u£ /or a radiatively induced first-order phase transition in a model with initial 
potential given by eq. fy.Q ) with (f>Q ka = 1.712 ko, X ko = 0.01 ko and g ko = 0.2 ko- The calculation is 



performed between the scales k{ 
units ofkf. 



-4.7 



ko and kf 



-5.2 



ko- All dimensionful quantities are given in 



In fig. 2 we consider a theory invariant under the symmetries <p\ <-> —<f>i, 4>2 ^ — 02> (f>i ^02- The 
potential is given by 



<^0 



u k0 (<pufo) = '-f [(4 - <g ko f + (4 - 4 k0 ?] + g -fm + v -^ [(4 - 4 ko f + (4 - 4J 3 } - (4-2) 

with 0Q fc() = 2 ko, Afc = 0.4 ko, g ko = 0.3 ko and u ko = 1. In fig. 2a we present the evolution 
of Uk((j>i) and in fig. 2b that of £/22(0i)> in complete analogy to figs, la and lb. The quantities 
W lk (r) = Un(Mr)) ~ U n (0) and W 2k (r) = U 22 (Mr)) ~ U 22 (0) are plotted in figs. 2c and 2d 
respectively. Our results for the nucleation rate are shown in fig. 2f. Again, we observe a significant k 
dependence of S k , ln(Aik / k'j) and ln(^2fc) which cancels out in the total rate. A qualitative difference 
between figs. 1 and 2 concerns the sign of ln(^42fc)> which is positive in the first case and negative 
in the second. This behaviour can be understood through the comparison of the form of W 2k (r) in 
figs. Id and 2d. In the first case, the mass term of the 4> 2 fluctuations in most of the interior of the 
critical bubble is smaller than that in the exterior. As a result, the operator — d 2 + U 22 (4>\>{r)) has 
lower eigenvalues than the operator — d 2 + U 22 (0) and the prefactor A 2k in eq. ([0]) tends to enhance 
the total rate. The opposite is true in the second case. The (ft 2 fluctuations become more massive 
in the interior of the bubble and they tend to be suppressed. As these fluctuations are an integral 
part of the critical bubble, the total rate must be suppressed as well. Notice that the form of W\ k (r) 
guarantees that the <f>\ fluctuations always tend to enhance the nucleation rate. 

In fig. 3 we present our results for a potential given by eq. fl4.2p with 0Q fco = 2 ko, \ ko = 0.35 ko, 
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Figure 5: ^4s in jig. [J, 6ti£ /or a radiatively induced first-order phase transition in a model with initial 
potential given by eq. with 4>^ kQ « 0.4999 ko, Xk = 0.01 ko and g ko = 0.1 ko- The calculation is 

performed between the scales ki = e~ 9 ' 2 ko and kf = e~ 9 ' 6 /co. All dimensionful quantities are given in 
units ofkj. 



9ko = O-^ an d Vk = 0-8- The behaviour is similar to that in fig. 2, but the first-order corrections 
ln(Aik/kj) and | ln(^42fc)| become now comparable to Sk- This signals the breakdown of the expan- 
sion around the semiclassical saddle-point. This is confirmed by the significant k dependence of the 
predicted total rate, which indicates that the higher-order corrections cannot be neglected anymore. 
The examples of figs. 1-3 demonstrate the reliability and consistency of our approach in theories with 
two fluctuating fields. The present work confirms that the picture of tunnelling we have developed is 



not relevant only for the one-field case studied in refs. [13, 25|. On the contrary, it has a wide range 
of applicability to many physical systems. 

We now turn to the discussion of radiatively induced first-order phase transitions. An example 
can be observed in a theory defined through the potential 



h 



Uk^uH) = '-f [(# - <P' 2 0k0 r + (02 - 4k f\ + ^14, (4-3) 

with 0Q feo = 1.712 ko, Afc = 0.01 ko and gu = 0.2 ko- Our results for this theory are presented in 
fig. 4. In fig. 4a we plot a large part of the evolution of Uk(<fti). The initial potential has only one 
minimum along the positive ^i-axis (and the equivalent ones under the the symmetries <f>i <-> —<f>i, 
(j>2 «-> — 4>2, 4>i ^2 ) and a maximum at the origin. In the sequence of potentials depicted by dotted 
lines we observe the appearance of a new minimum at the origin at some point in the evolution 
(at k/ko ~ e -4,4 ). This minimum is generated by the integration of fluctuations of the 02 field, 
whose mass depends on (j>\ through the last term in eq. ( f4.3|) (the Coleman- Weinberg mechanism). 
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Detailed descriptions of this phenomenon within the framework of the effective average action have 
been presented in refs. [14, 15, 23. 26, 36 for a variety of models. In fig. 4b it can be seen that the 
mass term of the fa field at the origin turns positive at the same value of k. This is a consequence of 
the fa <-> fa symmetry of the potential. We calculate the nucleation rate using the potentials of the 
last stages of the evolution. The solid lines correspond to ki/ko = e -4,7 , while the line with longest 
dashes corresponds to kf/ko = e -5,2 . In figs. 4b-4e we observe that the mass of the fa fluctuations in 
the interior of the critical bubble is much larger than the other mass scales of the problem, which are 
comparable to kf. This is a consequence of our choice of couplings g/X = 20. Such a large ratio of 
g/X is necessary for a strongly first-order phase transition to be radiatively induced. Unfortunately, 
this range of couplings also leads to large values for the fa mass and, as a result, to values of | ln(^4 2 fc)| 
that are comparable or larger than the saddle-point action Sk, even though ln(Aik/kf) remains small. 
As a result, the saddle-point approximation breaks down and the predicted nucleation rate I /kf is 
strongly /c-dependent. 

It may be possible to obtain a convergent expansion around the saddle point by considering models 
with smaller values of g. In fig. 5 we repeat the calculation of the nucleation rate for a theory described 
by the potential of eq. (4.3) with 0q, 0.4999 ko, Xk = 0.01 ko and gk = 0.1 k^. Due to the smaller 
value of the ratio g/X = 10, we observe a more weakly first-order phase transition. This can be 
checked by considering the value of kf that sets the scale for the dimensionful quantities in figs. 4 and 
5. This scale is smaller by a factor e -4 ' 4 ~ 0.012 in fig. 5 than in fig. 4. Fig. 5f indicates that the 
expansion around the saddle point is more problematic in this case. Not only | h^^fc)] is larger than 
the saddle-point action Sk, but the prefactor ln(Aik/kf), associated with the fluctuations of the fa 
field, becomes now comparable to Sk- This behaviour of A\k for weakly first-order phase transitions 



has already been observed in refs. 25]. We have not managed to find any region of the parameter 



space that leads to a convergent saddle-point expansion for the nucleation rate. The generic behaviour 
resembles either the one depicted in fig. 4 for strongly first-order phase transitions, or the one in fig. 5 
for weakly first-order phase transitions. The implications for cosmological phase transitions, such as 
the electroweak, are discussed in the following section. 



5 Discussion and conclusions 



In this paper we generalized to theories of two scalar fields the analysis of refs. [13, 25|] for the consistent 
description of first-order phase transitions in the one-field case. We established the reliability of the 
method by making consistency checks for a variety of models described by different potentials. The 



conclusions of refs. [13, 25] were verified in this more general context: Langer's theory of homogeneous 
nucleation is applicable as long as the expansion around the semiclassical saddle point that dominates 
tunnelling is convergent. The total bubble- nucleation rate is not determined simply through the 
exponential suppression by the action of the saddle point. The fluctuations around the saddle point 
are important and generate a pre-exponential factor which can be significant even in the cases for which 
the saddle-point expansion is convergent. In the example of fig. 1 the logarithm of the pre-exponential 
factor is about 30% of the saddle-point action. 

If the prefactor becomes comparable to the saddle-point action, the expansion ceases to converge. 
In our approach this is signalled by the appearance of a strong dependence of the predicted nucleation 
rate on the coarse-graining scale at which the free energy of the system is evaluated. We found two 
types of situations in which this happens: 

a) The fluctuations of the field fa, that serves as the order parameter for the phase transition 
and varies along the saddle point, become significant in the case of weakly first-order phase 
transitions. Typically these fluctuations enhance the total rate and can even compensate the 
exponential suppression. Several examples of this behaviour were given in refs. JITJ, SfJ. The 
conclusions were reconfirmed by the study of the prefactor A\k in fig. 5. The reason for this 
behaviour can be traced to the form of the differential operators in the prefactor associated with 
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the field <pi (see eqs. eqs. ( p.6| )-( |3~7c| )). This prefactor, before regularization, involves the ratio 
det'{-d 2 + mj + W lk {r))/det(-d 2 +mf), with mf = U n (0) and W lk (r) = U n (Mr)) ~ Un{0). 
The function W\k(r) always has a minimum away from r = (see figs. 1-5), where it takes 
negative values. As a result the lowest eigenvalues of the operator det'(— d 2 + m 2 + W\k(r)) are 
smaller than those of det(— d 2 + m 2 ). The elimination of the very large eigenvalues from the 
determinants through regularization does not affect this fact and the prefactor A\k is always 
larger than 1. Moreover, for weakly first-order phase transitions it becomes exponentially large 
because of the proliferation of low eigenvalues in det'(— d 2 + m 2 + W\k(r)). In physical terms, 
this implies the existence of a large class of field configurations of free energy comparable to that 
of the saddle-point. Despite the fact that they are not saddle points of the free energy (they are 
rather deformations of a saddle point) and are, therefore, unstable, they result in a dramatic 
increase of the nucleation rate. This picture is very similar to that of "subcritical bubbles" of 



ref. [37]. In ref. [38] the nucleation rate was computed by first calculating a corrected potential 
that incorporates the effect of such non-perturbative configurations. The pre-exponential fac- 
tor must be assumed to be of order 1 in this approach, as the effect of most deformations of 
the critical bubble has already been taken into account in the potential. In our approach the 
non-perturbative effects are incorporated through the prefactor. Both methods lead to similar 
conclusions for the enhancement of the total nucleation rate. 

b) The fluctuations of the field </>2, that is orthogonal to <^>i and can trigger a radiatively induced 
first-order phase transition, can also become important. In the radiatively induced scenario, 
the mass of 4>2 depends on the expectation value of 4>i through a mixing term ~ g4> 2 4>2 m the 
potential. If the coupling g is sufficiently large, a second minimum of the renormalized potential 
can be generated through the integration of the <p2 fluctuations, even if the bare potential has 
only one minimum. In all cases where this was achieved in the models we investigated (figs. 4 and 
5 and results not presented here), we found that the pre-exponential factor is large compared to 
the saddle-point action. This is not surprising. The radiative corrections to the potential and the 
pre-exponential factor have a very similar form of fluctuation determinants. When the radiative 
corrections are large enough to modify the bare potential and generate a new minimum, the 
pre-exponential factor should be expected to be important also. More precisely, the reason for 
this behaviour can be traced again to the form of the differential operators in the prefactor. The 
prefactor associated with the field 4>2 involves the ratio det(— d 2 + 77i| + W2k(r))/ det(— d 2 -l-m 2 ,), 
with m| = ^22(0) and W2k{r) = U22 {4>b{r)) — ^22(0). In units in which </>&(r) is of order 1, the 
function W2k( r ) takes very large positive values near r = (see figs. 4,5). This is a consequence 
of the large values of g that are required for the appearance of a new minimum in the potential. 
As a result, the lowest eigenvalues of the operator det(— d 2 + m| + M / 2fc( r )) are much larger than 
those of det(— d 2 + m 2 ). This induces a large suppression of the nucleation rate. In physical 
terms, this implies that the deformations of the critical bubble in the 4>2 direction cost excessive 
amounts of free energy. As these fluctuations are inherent to the system, the total nucleation 
rate is suppressed when they are taken into account properly. 

An interesting question concerns the possibility to ameliorate the situation for the radiatively 
induced first-order phase transitions by integrating out completely the <p2 field along the lines of ref. |8| . 
One may even envisage theories in which the term ~ §\ is absent and the gaussian integration over 
(f>2 can be peformed exactly. The result would be the appearance of a term ~ lndet(— d 2 + U22{4>i)) i n 
the "effective" action of the (f)\ field. For homogeneous field configurations, such a term would generate 
the two-minimum structure in the effective potential. However, the determination of the bubble 
profile would require the minimization of the whole "effective" action, including the additional term 
~ lndet(— d 2 + U22{4 ) i)) ■ The evaluation of the prefactor would involve a contribution that corresponds 
to a higher order in the semiclassical expansion that we employed in this work. The reason is that 
lndet(— d 2 + ^22(^1)) must be expanded around the saddle point configuration 4>b{r) and the new 
4>\ integration must be performed. Thus, it is clear that this approach does not offer any technical 
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advantages. Moreover, it is possible that an action that includes the term ~ lndet(— d 2 + ^22(^1)) is 
problematic due to the appearance of non-localities for field configurations for which U22(4'i) ~ 0. 

Our approach is based on a separation of the problem in two steps. At first, the high frequency 
modes of all fields are integrated out. The resulting effective theory can be approximated by a local 
action, as the integration of modes is performed with an explicit infrared cutoff so that non-localities 
are suppressed. Subsequently, the nucleation problem is addressed in a systematic expansion around 
the saddle point. A possibility suggested by the discussion in the previous paragraph is that an 
"improved" saddle point may be obtained if we minimize Sk + In Ak instead of only Sk- However, the 
corrections to the total rate that are induced by this change are of the same order as corrections from 
going beyond the gaussian approximation in the expansion around the saddle point. This has been 
confirmed by approximating InAk through eq. ( 3.14| ) and minimizing Sk -\-lnAk- For | ln^4fc| smaller 
than Sk we have found that, although the saddle point is modified, the value of Sk + ^nAk remains 
largely unaffected. 

Our results are relevant for most cosmological phase transitions. For example, the electroweak 
phase transition for Higgs boson masses below ~ 80 GeV is a radiatively induced first-order one. 
The same is true for the phase transitions in extensions of the standard model or in grand unified 
theories. Usually the second minimum of the potential of the scalar Higgs field appears as a result 
of the integration of the fluctuations of gauge fields, whose mass is generated through the Higgs 
mechanism. There is a close resemblance with the two-scalar models we considered. This is apparent 
if one considers the mass term of the 4>2 field along the 0i-axis. For the potential of eq. (|4.3| ), it is 
initially given by M| fco (0i) = -\k <f>ok / 2 + 9k 4>i/2- For the range of couplings g ko > A fco , which is 
relevant for radiatively induced first-order phase transitions, the first term in the previous expression 
is negligible. The evolution of the potential at scales k < ko results in the replacement of gk by an 
effective running coupling gk((fii). However, in cases similar to the ones we considered, in which no 
fixed points affect the evolution, the running of g is not substantial. 



In the gauged Higgs models, the mass of the gauge fields is generated through a term ~ e 2 cf) 2 A^A 



At the level of the bare action the mass term is M 2 ko (cpi) = ef, o (p 2 . Within a truncation that preserves 
the same invariants as the bare action, the evolution of the potential is given by an equation completely 
analogous to eq. (^^) with M 2 replaced by M 2 (and the modification of some numerical factors) [ |39| , 
p3|, |36|, 13(3] . The replacement of e? by a running coupling e|((^i) at scales k < k$ does not modify 



-k 

the qualitative behaviour of the potential, as long as a first-order phase transition takes place ]23|. 
It is very likely, therefore, that our conlcusions concerning the convergence of the expansion around 
the saddle point remain valid in the case of the radiatively induced first-order phase transitions in the 
gauged Higgs models. This would imply that the estimates of the nucleation rate that rely only on 
the exponential suppression factor may be very misleading. Moreover, the whole series may not be 
convergent and this would make an alternative calculational scheme necessary. 

A confirmation of the above possibility is provided by estimates of the bubble-nucleation rate based 
on a perturbative approximation of the effective action. Such a calculation was performed in ref. [II 
for the electroweak phase transition. An effective action for the Higgs field was obtained by integrating 
out the gauge fields through a perturbative calculation. The saddle-point action was evaluated within 
the effective scalar theory in the standard way, while the pre-exponential factor was computed using 
the heat-kernel method. There are several points that need refinement in this approach: a) It is 
questionable whether the gauge fields can be integrated out perturbatively, as they are massless in the 
symmetric phase within perturbation theory, b) The derivative expansion of the effective action for the 
scalar fields is problematic around the origin of the potential, where the gauge theory becomes strongly 
coupled, c) In the absence of a coarse-graining scale the fluctuations of the scalar fields cannot be 
treated properly, as the problem of double-counting their effect in the effective action and the prefactor 
is not resolved. It is remarkable, however, that, despite the above problems, the conclusions of ref. [|llj 
are in agreement with the implications of our results: The pre-exponential factor is significant and 
results in the additional suppression of the total rate. Moreover, the saddle-point approximation is 
problematic near the critical temperature. 
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The application of our formalism to the electroweak phase transition must deal with several techni- 
cal points. Firstly, an efficient truncation of the average effective action must be developed, appropriate 
for the discussion of gauged Higgs theories. Such truncations have been discussed in refs. |39], 23. 36, 40 



and the resulting evolution equations for the potential and the gauge coupling have been obtained. 
In the case of the SU(2) Higgs model, they account for the growth of the gauge coupling near the 
origin of the potential and the emergence of a strongly coupled regime. They also take into account 
the contribution of the Higgs-field fluctuations. The saddle-point action can be evaluated along the 



lines of refs. [14, 15, 13] and the present work. The calculation of the pre-exponential factor is more 
difficult. The fluctuation determinants associated with the gauge fields must be computed in com- 
plete analogy to those of the </>2 field in this work. This means that they must be evaluated in the 
background of the saddle point, with an ultraviolet regularization that matches the cutoff procedure 
in the derivation of the evolution equations for the potential and the gauge couplings. Particular care 
is required for the treatment of additional zero eigenvalues that result from the presence of Goldstone 
modes in the phase with symmetry breaking [41|. However, we believe that the above technical points 



should not affect the qualitative conclusion reached in this work: If the radiative corrections can affect 
the vacuum structure of the theory, the pre-exponential factor in the tunnelling rate is significant. 

Another possible application concerns the bound on the Higgs-boson mass from the requirement 
that our vacuum is not destabilized by the presence of a non-standard vacuum at large Higgs-field 
values . This new vacuum is generated by the top-quark radiative corrections to the effective 

potential. The most stringent bound arises from the requirement that the probability of thermal 
tunnelling from the symmetric phase to the non-standard vacuum at some temperature above the 
critical temperature of the electroweak phase transition is negligible. This problem can be treated 
within an effective three-dimensional theory at energy scales below the temperature. At such scales 
the top quark can be integrated out, as its propagator contains an effective mass of the order of the 
temperature. Its effect can be incorporated in the initial conditions for the evolution equations for 
the potential and the gauge couplings. The potentials that are relevant for this problem have a form 
very close to that along the spinodal line (the false vacuum is at the end of the metastability range 
and very close to becoming unstable) [43[|. It is, therefore, conceivable that the prefactor generated by 



the Higgs-field fluctuations is significant [25]. The effect of the gauge-field fluctuations must be taken 
into account as well, in analogy to the calculation of the nucleation rate for the electroweak phase 
transition. 
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